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Two applications of mean-field calculations based on 3D coordinate-space techniques are presented. The first 
concerns the structure of odd- TV superheavy elements that have been recently observed experimentally and shows 
the ability of the method to describe, in a self-consistent way, very heavy odd-mass nuclei. Our results are 
consistent with the experimental data. The second application concerns the introduction of correlations beyond 
a mean-field approach by means of projection techniques and configuration mixing. Results for Mg isotopes 
demonstrate that the restoration of rotational symmetry plays a crucial role in the description of 32 Mg. 



In this note, we shall first present a very representative 
application of the method to odd-N superheavy elements 9 ) 
(SHE). We shall then describe new developments 10 -' re- 
lated to the restoration of the rotational and particle 
number symmetries; these improvements are expected to 



Introduction 

Mean-field calculations on a 3-dimensional mesh were 
introduced in nuclear physics in the late 1970s in the 
studies of nuclear collisions within the time-dependent 
Hartree-Fock method 1 ) . In these early applications, schematise extremely useful in future studies of medium-mass nu- 
nuclear interactions were used, without spin-orbit and clei far from stability. 
Coulomb terms, and applications were limited to colli- 
sions between light nuclei. The same technique has later 
been applied in a systematic study of the spectrocopic 
properties of Zr and neighboring isotopes 2 ). New de- 
velopments have then permitted the use of effective in- 
teractions of the Skyrme type, including Coulomb and 
spin-orbit terms, and pairing correlations were considered 
the BCS approximation with a seniority force. There 
are several advantages of the coordinate-space technique 
which make it very attractive. It allows one to describe 
very general intrinsic deformation with a good numerical 
accuracy. It also provides a very flexible computation 
scheme, suitable for generalizations (e.g., introduction of 
correlations beyond the mean-field theory). The main 
developments of the method performed after the comple- 
tion of the Zr study were the introduction of a cranking 
constraint to describe rotating nuclei 3 ) , the introduction 
of the configuration mixing using the generator coordi- 
nate method 4 - 1 , the particle-number projected HF+BCS 
calculations 5 ' , and improvements in the treatment of the 
pairing correlations by the use of the HFB method 6 ) 
with density-dependent zero-range interactions 7 ). Let 
us also mention a theoretical study of the properties of 
coordinate-space calculations 8 ) which introduced power- 
ful ways to interpolate wave functions discretized on a 
3D mesh. 



Structure of Odd-iV Superheavy Elements 

There is no consensus among theorists with regard to 
the center of the shell stability in the region of spher- 
ical SHE. For the neutrons, most calculations predict 
a magic gap at 7V=184. However, because of different 
treatments of the large Coulomb potential and spin-orbit 
interaction, various models yield different predictions for 
the position of the magic proton gap (cf. discussion in 
Ref. 11 '). In this context, the synthesis of the iV=175 
isotones reported by the Dubna/Livermore (DL) collab- 
oration (Z=114) 12 ) and the Berkeley/Oregon (BO) team 
(Z=118) 13 ) is particularly important. Due to the large 
neutron number of the compound nucleus, the observed 
a chains do not terminate at some known systems. Con- 
sequently, until the charge-mass identification is made, 
one can use theoretical arguments to support/disprove 
the experimental assignment. 



We have performed a large-scale self-consistent study of 
properties of even-even, ^4-odd, and odd-odd heaviest 
nuclei. In our work, we have used the Hartree-Fock- 
Bogoliubov (HFB) method with a Skyrme interaction in 
the particle-hole (ph) channel and a delta force in the 
pairing channel. The details of the calculations closely 
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Table 1. Predicted structure of the lowest one-quasi- 
particle excitations in the 289 114 a-decay chain. Each 
excitation is characterized by the Nilsson quantum 
numbers of the odd neutron, excitation energy, and 
quadrupolc deformation. The number in parenthe- 
ses indicates the Q a value for the ground-state to 
ground-state transition. 
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follow Ref. n \ The HFB equations have been solved in 
the coordinate space according to the method of Ref. 7 - ) . 
In the ph channel, the Skyrme effective interaction SLy4 
has been used. This parameterization 14 ) has, in particu- 
lar, been adjusted to reproduce long isotopic sequences; 
hence, one can expect it to have good isospin proper- 
ties. The pairing strengths of the interaction have then 
been adjusted to reproduce the average proton and neu- 
tron pairing gaps in even-even nuclei around 2 Fm. The 
properties of one-quasiparticle states were calculated by 
means of the self-consistent blocking. 

The calculated one-quasiparticle structures and the g.s.- 
to-g.s. values of Q a in the 289 114 and 293 118 a-decay 
chains are shown in Tables 1 and 2, respectively. 

According to our calculations, the g.s-to-g.s. a decays of 
289 114 and 293 118 are structurally forbidden since the g.s. 
properties of parent and daughter nuclei differ dramat- 
ically. The allowed transitions to the excited [707] 4^ 
level in 285 112 and 289 116 correspond to Q Q =9.0 MeV 
and <2a = 10.9 MeV, respectively, and are considerably 
lower than the experimental values (DL: 9.9 MeV; BO: 
12.6 MeV) Consequently, the most probable candidates 
for the first a transitions are the 10.2 MeV ( 289 114) 
and 11.8 MeV ( 293 118) lines associated with the allowed 
[611]^ + ^[611]^ + decays. By the same token, the [611]| + 
g.s. of V=173 isotones is expected to decay to the ex- 



Table 2. Same as in Table 1, except for the 293 118 



a-decay chain. 
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cited [611]i level in the V=171 daughters. For 285 112, 
the corresponding Q a energy, 8.76 MeV, is very close to 
the experimental energy of the second a transition, 8.84 
MeV, reported by the DL group. For the 289 116 decay 
we obtain Q Q =10.14 MeV; i.e., we underestimate the en- 
ergy of the second a particle (11.8 MeV) from the BO 
experiment. This is the worst deviation from experiment 
obtained in our calculations. 

The [611] 2" + level is also expected to be the ground state 
of the 7V=169 isotones, and our prediction for the al- 
lowed [611]i + ^[611]i + decays is 9.4 MeV for 281 110 
(DL: 9.0 MeV) and 10.6 MeV for 285 114 (BO: 11.5 MeV). 

For the a transitions in the 281 112^ > 265 104 chain, 

allowing the positive-parity Nilsson orbitals with very 
similar quantum numbers, we obtain the following val- 
ues of Q a : 10.7 MeV (BO: 10.8 MeV), 11.1 MeV (BO: 
10.3 MeV), 9.9 MeV (BO: 9.9 MeV), 8.0-8.5 MeV (BO: 
8.9 MeV). An alternative route is possible that involves 
a transitions between j'15/2 orbitals. Here, for the last 
two Q a values we obtain 9.8 MeV and 8.7 MeV. In both 
cases we obtain good agreement with BO data. 

The calculated equilibrium deformations of one-quasiparticle 
states in V=175 isotones are rather small (/?2~0.11), and 
they increase along the a-decay chain. This trend re- 
flects the influence of the V=184 magic neutron gap for 
the heaviest systems and the deformed V=162 gap for 
the lightest nuclei in the chain. As far as Q a values are 
concerned, our calculations are consistent with the recent 
experimental findings. 
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Correlations beyond mean field in Mg isotopes 

The cranking method is widely used in nuclear spec- 
troscopy to describe high-spin states. Applications based 
on effective nuclear interactions have been particularly 
successful in the description of superdeformed rotational 
bands in several regions of the mass table. In the crank- 
ing method, a rotational band is generated by the rota- 
tion of a deformed intrinsic state. Since cranking states 
are not eigenstates of the angular momentum, this causes 
some problems in determining, e.g., transition rates in 
nuclei which are not very well deformed. 

Another limitation of the cranking model appears in nu- 
clei which are soft with respect to the variation of a collec- 
tive variable. In this case, one expects the interference 
of the zero-point vibrational mode with the rotational 
motion which leads to variations in the nuclear structure 
along the yrast line. 

In ref 5 ) , we have presented a method to restore the par- 
ticle number symmetry within the HF+BCS theory that 
allows us to perform a configuration mixing of projected 
wave functions in the direction of selected collective vari- 
ables. The method presented in this section generalizes it 
by the inclusion of a restoration of the rotational symme- 
try. It enables us to describe details of collective spectra 
and transition rates. 

The starting point of the method are wave functions I^q) 
generated by mean-field calculations with a constraint 
on a collective coordinate a. Wave functions with good 
angular momentum and particle numbers, 
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are obtained by means of the projection operators P. 
In Eq. (|l|), N x is a normalization factor, depending on 
x = J, M, N, Z, a. 

Using the projected state |<I>, JMa) as a generating func- 
tion, configuration mixing along the collective variable a 
is performed for each angular momentum: 



|*,JM)=^/ Q 7M |$,JMa). 



(2) 



The weight functions /„ M are found by requiring that 
the expectation value of the energy, 



E 



JM = (%JM\H\%JM) 
(*,JM|*,JM) ' 
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This prescription leads to the discretized Hill- Wheeler 
equation 
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in which the Hamiltonian kernel H JM 
kernel T JM are defined as 
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and the overlap 
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The kernels are obtained by integrating the matrix el- 
ements between rotated wave functions over three Eulcr 
angles and two gauge angles. Besides these kernels, one 
can determine transition probabilities between different 
eigenstates of the Hill- Wheeler equation. This requires 
the calculation of matrix elements of electromagnetic op- 
erators. 

Currently, in order to save computing time and to test 
the method, certain symmetry restrictions have been im- 
posed on the mean-field wave functions. Namely, they 
have been assumed to be axially symmetric and time- 
reversal invariant. In this way, the integration over the 
Euler angles is limited to a single angle. Pairing correla- 
tions are treated in the BCS approximation. 

In the following examples, a Skyrme force is used in the 
particle-hole channel and a density-dependent zero-range 
interaction in the pairing channel. As in ref 5 ), in the 
calculation of non-diagonal matrix elements, the density 
dependence of the Skyrme interaction is generalized to a 
dependence on the mixed density. 

Application to 2A Mg 

The results shown in this section have been obtained us- 
ing the HF+BCS wave functions generated with an axial 
quadrupole constraint. The Lipkin-Nogami prescription 
has been used to improve the treatment of pairing cor- 
relations. The variation of the energy as a function of 
prolate and oblate deformations is plotted in Fig. 1 for 
the Sly4 Skyrme interaction and a surface pairing inter- 
action having strength G=1000MeV fm 3 for both neu- 
trons and protons. In addition to the prolate absolute 
minimum, the mean-field curve presents a shoulder at an 
oblate deformation around 50 fm 2 . The triple projection 
creates an oblate minimum at the position of the shoul- 
der for J up to 6 + . For greater values of J, the weights 
of the intrinsic wave functions for deformations below 
-200 fm 2 are very small. Consequently, the projected en- 
ergy curves do not exhibit any oblate mimima. 

For each value of the angular momentum, we have per- 
formed a configuration-mixing calculation including quadrupole 
moments between -350fm 2 and 450 fm 2 . This corresponds 
to intrinsic configurations excited by about 30 MeV with 
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respect to the prolate minimum. The spectrum that is 
generated in this way (represented by bars) is plotted 
at the quadrupole moment corresponding to the largest 
component of the collective wave function. The value of 
this quadrupole moment is very close to the minimum of 
the projected energy curve. Moreover, the energy of this 
minimum is only slightly modified by the configuration 
mixing. The largest gain, ~800keV, is obtained for the 
+ state, but it is reduced at higher spins. Several ex- 
cited states are found at low energy for each spin value. 
Except for the second + and 10 + , the wave functions 
of yrare states are peaked around the oblate secondary 
minimum. 

In the same figure, theoretical transition probabilities 
along the yrast line are compared with the data. In 
the GCM calculation, the minima of projected energy 
curves have been used. The B(E2; 2f — ► Of) rate cal- 
culated between the minima of the two curves is very 
close to the experimental value. The configuration mix- 
ing causes a spreading of the collective wave function 
on the quadrupole moment and decreases the value of 
the B(E2) rate. Since for spins greater than zero wave 
functions do not have low quadrupole-momcnt compo- 
nents, the configuration mixing docs not affect signifi- 
cantly the transition probabilities. Here, the agreement 
between calculations and experiment is excellent. One 
has to note, however, that the inclusion of triaxiality is 
expected to modify these predictions to some extent. 

Application to 32 Mg 

Altough Mg corresponds to the A^=20 shell closure, 
there is some experimental evidence that it is deformed 
in its ground state. The excitation energy of its first 2 + 
state is only around 900 keV, i.e., significantly lower than 
in lighter Mg isotopes. Also the B(E2;2f — ► 0+) value 
is unusually large. 

Figure 2 shows the calculated energy curves for 32 Mg. 
The ground state calculated in the HF+BCS approxima- 
tion is spherical (see discussion in ref. 15 ). However, unlike 
in the 24 Mg case, the restoration of broken symmetries 
modifies the topology of the energy curve in a more dra- 
matic way. Namely, the resulting + ground state cor- 
responds to the projection of a deformed intrinsic state 
with a quadrupole moment of 100 fm 2 . The minima for 
the higher- J values can be associated with even larger de- 
formations. The configuration mixing leads to small en- 
ergy gains. The collective + wave function is spread over 
a large range of quadrupole moments between -200 fm 2 
and 200 fm 2 , while for other J- values the weights of the 
HF+BCS wave functions have well-pronounced maxima 
around 200 fm 2 . 

The calculated spectrum and transition probabilities are 
compared to the experimental data in Fig. 2. It is seen 
that our calculations do not reproduce the pronounced 
collectivity of 32 Mg. Compared to the data, the energy 



of the 2+ is too high and the B(E2) value is too low. It 
should be noted, however, that the delicate balance be- 
tween the spherical minimum and the deformed intruder 
configuration in this nucleus is sensitive to the effective 
interaction used 15 ). 
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Fig. 1. Symmetry- projected results for the nucleus 24 Mg. (a) Energies as a function of the axial quadrupole 
moment. Solid line with stars: HF+BCS+LN results; the curves displaying filled circles, filled boxes, filled 
diamonds, circles, boxes, and diamonds correspond, respectively, to projected energies in the J=0 + , 2 + , 4 + , 6 + , 
8 + , 10+ states. The horizontal bars are the projected energies of the collective intrinsic states, (b) Excitation 
energies and transition probabilities (numbers in e 2 fm 4 ) between: intrinsic states minimizing each projected 
energy curve (left), yrast collective states (right), and experimental levels (center). 
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